Dynamical symmetries and well-localized 
hydrogenic wave packets * 

Vladimir ZVEREV f and Boris RUBINSTEIN * 

' Ural State Technical University, Ekaterinburg, 620002, Russia 
zverev@dpt.ustu.ru 

* Department of Mathematics, University of California, Davis, CA 95616, USA 

boris @math. ucdavis. edu 

Abstract 

A new method for constructing of composite coherent states of the hydrogen atom, based 
on the dynamical group approach and various schemes of reduction to subgroups, is presented. 
A wide class of well-localized (Gaussian) hydrogenic wave packets for circular and elliptic orbits 
is found using the saddle-point method. 

1 Introduction 

In recent years, new experimental techniques opened way to creation and study of high energy 
(Rydberg) states in atoms. These states are described by approximate hydrogenic wave functions 
with very large principal quantum numbers. Some new effects, as the dynamical localization and 
the dynamical chaos, have attracted considerable interest. Explanation of these phenomena uses 
classical equations of motion [T] . It is reasonable to look for an alternative quantum description on 
the basis of semi-classical approximations, which is naturally provided by a coherent states (CS) 
formalism. 

In Section 2, starting from the 0(4,2) dynamical group approach [2] and using three schemes 
of reduction to subgroups 0: 0(4,2) D 0(4) - 0(3) ® 0(3), 0(4, 2) D 0(2,2) - 0(2,1) <g> 
0(2, 1), 0(4, 2) D 0(3) ® 0(2, 1), we construct composite CS in physical and auxiliary ("tilted") 
representations 0). We use two types of generating operators of CS with different procedures of 
transition to a classical limit. In particular, the generating operators for Perelomov S0(3) and 
S0(2, 1) CS 0, Barut-Girardello S0(2, 1) CS HE], generalized hypergeometric CS |7j, Brif S0(3) 
and S0(2, 1) algebra eigenstates [8], may be used for this purpose. The CS are separated into two 
classes with different semi-classical behavior. 

The hydrogenic CS wave functions have complicated form, so it is reasonable to use simplified 
asymptotic expressions. In Section 3 we describe a method for asymptotic estimate and obtain well- 
localized hydrogenic wave packets for circular and elliptic orbits. A similar asymptotic estimate 
method is used in the theory of CS path integrals. We believe that the approach discussed in this 
paper can be applied to computation of the CS path integrals for the hydrogen atom and other 
systems with known dynamical symmetry. 



2 The hydrogen atom: 0(4, 2) dynamical group and reduc- 
tions to subgroups 

Denote generators of the dynamical group 0(4, 2) of the hydrogen atom 00] (the group of rotations 
in six-dimensional pseudo-Euclidean space with a metric g = diag(l, 1,1,1,— 1,-1)) as Cap (in the 
notation used in |3]): 

£ a0 = (r x p) a , C al = (r a p 2 - 2p Q rp + r a )/2, 
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£ Q 2 = -rp a , C a3 = (r a p 2 - 2p a rp - r Q )/2, (1) 
£01 - (rp 2 - r)/2, £ 02 = rp - i, £ 03 - (t 2 + r)/2. 

The Schrodinger equation for the hydrogen atom reads (% — = with Ti = p 2 /2 — a/r . 

After multiplication by r it can be rewritten in terms of the generators Q ED E3 EOj : 

r(W - E)\V) = {{C 03 + £oi)/2 - S(£ 03 - £ i) - o}|*) = 0. (2) 

The unitary transformation 

I*) = cxp(?;0£ O2 )|*), tanhfl = (1 + 2£)/(l - 2E), 

reduces to 

(£ 03 -a/\/=2£)|¥) =0. 



The eigenstates of £03 may be chosen in this auxiliary representation (AR) as \n, I, m) sp h or 
\n, ni, n 2 , m) par with n = a/y/—2E, in the spherical or parabolic coordinates. An inverse transfor- 
mation 

\n,l,m)s P h l- CX p H 0£ O2 )/ \n,lM»h 



n,ni,n2,m)p ar J [ \n,n ll n2,m) par 

produces the physical representation (PR) of the eigenvectors. 

Consider three types of reduction to subgroups and corresponding Lie algebras: 

(a) 0(4,2) D 0(4) -0(3)® 0(3), o(3) © o(3) = {V^} © {7><f>}, 

(6) 0(4, 2) D 0(2, 2) ~ 0(2,1)® 0(2,1), o(2, 1) © o(2, 1) = {Q^} © {Q^}, (3) 

(c) 0(4, 2) D 0(3) 0(2,1), o(3) © o(2, 1) = {£ a0 } © {£ 0q }, 

where 

= (£aO ± C a3 )/2, Q&) = (£ a ± C 3a )/2. 
For these cases the space R of bound states of the hydrogen atom is decomposed as follows: 

(a) R=® n > 1 {R(°Wn-l,(n 1 ))®R( Wn-l,(n 2 ))}, (4) 
(6) R = ® meZ {R{°Wl + \m\, (ni + ^(m - |m|))) ® + Hj (n2 + * (m _ H)))}; 

(c) i? = ®,> {i?(°( 3 >2Z, (Z - m)) © R{° {21) 21 + 2, (n - 1 - 0)}, 

where R( G S, (K)) — { | G 5, (if)) } denotes the irreducible representation space of a group G = 0(3) 
or 0(2, 1) with generators {E a }, for which we have 

E 3 \ G S, (K)) = (S/2-e G K)\ G S, (K)), 

{El + Ei + e G E 2 )\ G S, (K)) = (S/2 - e G S 2 /A)\ G S, (K)), 

where £o(3) — 1) e o(2,i) = ~L an d the representation space R( G S, (K)) is spanned by the set of 
orthonormal eigenvectors \ G S, (K)) . The vector subspaces 

R [ +] = {\n,n- l,ra- l) sph } = {\n, 0,0, n- l) por }, 

i?^ 1 = {|n,n- 1, -n + l) sp fc} = {|n,n- l,n- l,-n + l) por }, 

with maximal m = n - 1 or minimal m = — n + 1 values of the magnetic quantum number 
correspond to the circular orbits The vectors from R^ play the role of reference vectors in 

decompositions (@J). On the other hand, R Q ' and i?Q + ' are irreducible spaces of 0(2, 1) (S = 1). 
The corresponding 0(2, 1) Lie algebras have the form 

nf ] = (±£n - £22)72, n l 2 ±] = (£21 ± £12)72, n{ ±] = (±£ 03 ± £ 30 )/2, 
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where ± corresponds to the subspace R with m — ±(n — 1). Following 12, we construct 

coherent states in two stages. At hrst, using the operators {7i|^}, we construct the circular orbits 

CS in Next, we construct CS in the global space R using one of the generator sets from 

Each reduction scheme corresponds to a specific construction method of the hydrogen atom CS: 

(a) |°( 4 )^M,r, ?? )=I? I ( ( 3 )r,(^ + )))P I ( ( 3 )r;,(^-)))|c[ ± l), 

(b) \°™ u Wt, V ) = Th^r, (Q^)W°^ V , {QtV% (5) 

(c) ^ffl^P.^W.^l^.rr,^))^^^^))!^^), 

where 

\J ±] ) =Pn( o(2a ^,(HL ± l))|l,0,0,0) par (6) 
is the circular orbit CS. We use an uniform notation for the generating operator of CS: 

\ G S, (u)) = Vn( G uj,(E a ))\ G S, (0)), 

where G — 0(3) or 0(2, 1), and f2 = I or II. An important point is that one should use two 
types of generating operators of CS, with different procedures for transition to the classical limit: 
S -> oo (n = I); M -> oo (n = II). 

There exist many algorithms to construct coherent states of dynamical groups and quantum 
physical systems (see review ^3]). In most cases these CS provide a transition to the "classical 
limit" in some sense. However, a selection of the CS for a specific physical system must be done 
with great care. It is essential to take into account not only the symmetry of a model, but also 
the type of semiclassical behavior. We consider some special types of CS and clarify conditions of 
quantum-classical correspondence. 

Usually CS |£) is called semiclassical if an overlapping distribution function |(£ |£ + S^)\ 2 has a 
sharp peak for small |<5£|, and becomes singular delta-shaped function in the classical limit, which 
can be written in the form 

<£ |„) cl ^- o,t?m <£ 10 = i. (?) 

Satisfaction of the condition 

(Z\E 2 a \0/mE a \0) 2cl ^l, (8) 

where E a are elements of the Lie algebra L, gives another evidence of semiclassical properties of 
the CS The conditions 0, (jHJ) are satisfied in the cases: the Barut-Girardello 0(2,1) CS, 
|£| -► oo, L = o(2,l); the Perelomov 0(2,1) CS, S ^ oo, L = o(2,l); the Perelomov 0(3) CS, 
S — ► oo, L = o(3). 

An extended class of semiclassical CS follows from a definition of the generalized hyper geometric 
CS 0: 

where (a) n is the Pochhammer symbol and p F q denotes the generalized hypergeometric function; 
£ is a complex argument; integers Q!j are negative for i = I,..., I and positive for i = I + 1, ...,p; 
< I < p; pj are positive real numbers. The overlapping function for this states has the form: 

<£|C) - N^N^ p F q ((a p ), { Pq ), {-l) l C0, m = I- 

Starting from the definition @, let us consider a number of special realizations: 

(i) p = q + 1 > 1, pj = 1, a.j — S, S — 1, 2, ...; in the special case p = 1 the state (0 is the 
Perelomov 0(2, 1) CS 0; 

(ii) p = q + 1 > 1, pj = 1, Oj = —5, 5 = 1,2, ...; in the special case p = 1 the state @ is the 
Perelomov 0(3) CS 0; 

(iii) q > p, pj — 1, otj = S, S = 1, 2, ...; in the case p = q = the state @ is the CS of a 
harmonic oscillator |14| : in the case p = — 1, q = this is the Barut-Girardello CS 6 . 
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It can be shown directly that for the cases (i) and (ii) the requirement Q is satisfied for the 
limiting condition S — > oo (ft = I), while for the case (iii) the proper condition is |£| — > oo (Jl = II). 

For comparison, we consider briefly an alternative extension of a family of CS for 0(3) and 
0(2,1) groups - the Brif's algebra eigenstates 8: 

+ P 2 E 2 + 3 E 3 ) \C,Pu P2,fe) = C |C, Pi , ft, fa) , 

where |£, /3i, ft, A) is a linear superposition of \ G S, (if)) with different if. In the general case 
6 = + /3| + e G /?I) 1/2 ^ and for C = (/ - e G S/2)b, 1 = 0,1, 2--, the condition Q is satisfied 
at S — > oo (O = I). At the same time, in the degenerate case 6 = the 0(2, 1) algebra eigenstates 
are semiclassical at the both S* oo and |C| — + oo limiting conditions (0 = 1 and II). 
These examples do not cover all diversity of CS with semiclassical behavior. 



3 Semiclassical asymptotics of coherent states wave func- 
tions 

We use the standard saddle-point method of asymptotic estimate for integrals 

F(A,x) =^exp{A/(x,t)}0(x,t)dt, 

where A > 1, x = {x±, x%, ■ ■ ■ , x n } and t = {t\, t%, i m }; /(x,t) and </>(x,t) are holomorphic 
functions; c is a n-dimcnsional smooth manifold deformed to reach a minimax (a saddle point is 
nonsingular). The asymptotic formula reads |15j 

F(A,x) ~ (2^/A) m/2 exp{A/(x, t o (x))}0(x, t (x))/V dot [—fi (x, t (x))], (10) 

where t = to(x) is a solution of the system of equations df/dti = 0, i = 1, m, for a stationary 
point, and /Lt«j = d 2 f /dUdtj. We also assume that |F(A, x)| has the delta-shaped peak near the 
point x = x o- Since |F(A, x)| ~ | exp{A/(x, t (x))}| ~ exp{ARe/(x, t (x))}, we can find x o as 
an extremal point of the real part of the exponent. The supplementary equations for the stationary 
point have the form 

a/(x,t (x)) , a/(x,t) | 

— ^ — + c - c - = -^r l t=t °« + cx - = °- 

Substituting the second order expansion of /(x, t (x)) at x = x o into IjlOfl . we find the desired 
expression for the asymptotic estimate 

F(A,x) ~ (2tt/ A) m/2 y/ det P/Xoo] Voo exp {A/ o + iAh T <5x 
+ (A/2) 6x T (ctoo ~ 7^0 i'oo 1 7oo) <5x} , 

with 



on 



where we use the notation (-) 00 = (-)\ x _ x t=t for different functions given in the point (xoo, too), 
too = to(xoo). 

Consider in detail a special case of hydrogenic CS wave function. Using the formulae J5J:) and 
10 with generating operators for the Perelomov CS (P) and the Barut-Girardello CS (BG), we 
obtain 

k 0, V ) = 2M£ (2,1) ?7, {Co a ))Vu{ B { G' X) u, (HL ±] ))|1, 0, 0) sph . 
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Noting that variation of r in implies only a rotation of the CS wave function without change 
of its shape, we can set r = 0. Using the coordinate hydrogenic wave function in the auxiliary 
representation 



(n + Q! 
(21 + 1)! V (n-l-l)V 



(r\n,l,m) = " J , v " 7 ' J ' (2rf e - r 1 F 1 {l + 1 - n, 21 + 1, 2r)y im ( 



and in the physical representation 



1 r 

r|n, Z, m) = — (— |n, Z,m) 



n 

we obtain the explicit form for the CS: 



Performing the summation we arrive at 

r \u},0,i]) = -2 exp -r- I , 

y/irl (2\v\) (1 - V) 2 \ 1-VJ V 1 ~V J 

where r + = rexp(ifa), r = |r| and Iq denotes the modified Bessel function. The asymptotic 
estimate (obtained directly or by the saddle-point method) for the absolute value of the CS wave 
function has the form: 



2^/ 2 \w\ 1 / 2 |l-rj| 4 *\ 2\u\ \1 - r)\ 
where Ax = x — xq, Ay — y — yo; the point with coordinates 



xq = (n)oo(e - cosd), y = (n)ooVi - e 2 sin( 



lies on the elliptic orbit with an eccentricity e: 



H(i + H 2 ) _ 2H 



i-M 2 ' i + H 2 ' 

In a similar way one can find the asymptotic estimate for the physical CS wave function, but it is 
omitted here due to its complexity. In this case only the saddle-point method with interchange of 
summation and integration can be used. It is remarkable that in this case the system of equations 
for the saddle point admits analytical solution. In a special case of circular orbits the asymptotic 
estimate for auxiliary and physical representations simplifies to 

* , v 1 f , , * , z 2 \lo\A6 2 Ap 2 iApAcb] , . 

(r ° + Ark °' 0> * ^/4 M i/4 exp {'H ^ ~ 2R " 4 - 4R + -IT } ' (n) 

where r = (pcos fa psin fa z), A(j) = <f> — 4>o, Ap = p — \uj\, and 

r + Ar u. 0, w -= — ; exp <^ i\uj\A6 - — — - J— L_T £_ + £l_r L 12 

V 1 ' ' 1 v^ 3 / 4 M 9 / 4 I 2|w| 3 10 10M 3 5|w| J y ' 

where r = (pcos fapsmfa z), A<f> = <j>— <f>o, Ap = p — \oj\ 2 . The functions and 1)12(1 are related 
by the following equation: 

(r 1 |w,0,0)= f ^■K(T 1 ,r 2 )(T 2 \ZJftfy, 

where 

1 /. z 2 2 z 2 p 2 A0 2 Ap 2 2zApA^\ 

and A(j) = fa - fa, Ap = p x - p\. 
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